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[i]  Combined  Release  and  Radiation  Effects  Satellite  (CRRES)  Electric  Field  Instrument 
(EFI)  data  are  used  to  determine  the  electric  field  power  spectral  density  as  a  function 
of  L  and  Kp  over  the  frequency  range  0.2  to  15.9  mHz.  The  power  at  each  frequency  is  fit 
to  the  function  P(L,  Kp)  =  a  Lb  exp(cKp).  Assuming  a  purely  electrostatic  field  and 
making  several  other  assumptions  regarding  the  azimuthal  dependence  of  the  field 
fluctuations,  a  Kp-dependent  radial  diffusion  coefficient  Dll  is  computed  from  the  power 
spectra.  The  model  average  Dll  for  high  activity  (Kp  =  6)  are  between  1  to  2  orders 
of  magnitude  larger  than  that  for  low  activity  (Kp  =  1),  dependent  upon  L  and  first 
invariant. 
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1.  Introduction 

[2]  Delineating  the  specific  acceleration  mechanisms  that 
populate  and  maintain  the  relativistic  electron  radiation  belt 
is  of  key  interest  to  the  magnetospheric  and  space  weather 
community.  There  has  been  much  attention  given  to  various 
acceleration  mechanisms  in  which  the  electrons  interact 
over  a  range  of  time  scales  with  either  MHD  pulses,  ULF 
waves,  or  VLF  waves  (for  a  review,  see  Friedel  et  al. 
[2002]).  Radial  diffusion  as  formulated  in  the  early  1960s 
[e.g.,  Falthammar ,  1965]  has  long  been  considered  a 
fundamental  acceleration  mechanism  that  gradually  (over 
several  months  time  scale)  energizes  the  radiation  belt 
particle  population.  Through  random  perturbations,  a  par¬ 
ticle’s  third  adiabatic  invariant  (related  to  its  azimuthal  drift) 
may  be  violated  while  its  first  two  invariants  (related  to  its 
gyro  and  bounce  motions,  respectively)  are  conserved. 
Through  such  a  stochastic  process,  electrons  gain  energy 
as  they  diffuse  towards  lower  phase  space  density  at  lower  L 
and  greater  magnetic  field.  A  dynamic  balance  of  radial 
diffusion  (source)  and  pitch  angle  diffusion  (loss)  leads  to 
the  observed  double  peaked  radiation  belt  radial  profile 
[Lyons  and  Thorne ,  1973].  More  recent  efforts  have 
extended  the  early  theoretical  work  on  radial  diffusion 
to  include  drift-resonant  interactions  between  electrons 
and  ULF  toroidal  [Elkington  et  al .,  1999]  and  poloidal 
[Elkington  et  al .,  2003]  wave  modes.  These  studies  have 
found  that  these  drift-resonant  interactions  can  lead  to 
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diffusion  rates  at  least  as  significant  as  those  predicted  by 
the  earlier  treatments. 

[3]  Estimates  of  the  radial  diffusion  coefficient  due  to 
electric  and  magnetic  fluctuations  (DLl  =  Dll  +  Dll)  have 
been  made  from  particle  observations  over  long  time 
intervals  (months)  [Frank,  1965;  Newkirk  and  Walt ,  1968; 
Lanzerotti  et  al .,  1970;  Lyons  and  Williams ,  1975;  West  et 
al .,  1981;  Selesnick  et  al .,  1997].  These  studies  presume  a 
single  DLl  fixed  for  all  times  independent  of  magnetic 
activity.  Other  studies  determine  Dll  and  Dll  fr°m  direct 
field  measurements.  Holzworth  and  Mozer  [  1 979]  analyzed 
electric  field  data  from  a  balloon  campaign  at  L  =  6  and 
found  a  Kp-dependent  Dll-  Mozer  [1971]  had  earlier  found 
that  the  power  spectra  related  to  Dll  show  a  Kp  dependence 
but  no  L  or  local  time  dependence.  Magnetic  field  power 
spectral  densities  and  the  corresponding  Dll  have  been 
determined  at  L  =  4  [Lanzerotti  and  Robbins ,  1973; 
Lanzerotti  and  Morgan ,  1973]  and  at  geosynchronous 
orbit  [Arthur  et  al .,  1978;  Lanzerotti  et  al .,  1978]. 

[4]  It  has  been  shown  that  an  activity  dependent  DLl 
spanning  orders  of  magnitude  is  necessary  (but  not  suffi¬ 
cient)  to  model  magnetic  storm  time  behavior  of  relativistic 
electrons  [Brautigam  and  Albert ,  2000].  However,  until 
now,  there  has  been  no  direct  evaluation  of  electric  field 
power  spectral  densities  throughout  the  equatorial  inner 
magnetosphere  (L  =  3.0-7.0)  for  a  wide  range  of  magnetic 
activity  (Kp  =  0  to  7)  that  could  be  used  to  directly 
determine  DlL(L,  Kp). 

[5]  The  Combined  Release  and  Radiation  Effects  Satellite 
(CRRES)  electric  field  instrument  (EFI)  [Wygant  et  al ., 
1992]  has  provided  the  opportunity  to  greatly  extend  our 
understanding  of  the  inner  magnetospheric  electric  field. 
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Previous  work  using  CRRES  EFI  data  has  addressed  the 
large  scale  DC  electric  field  [Rowland  and  Wygant ,  1998; 
Wygant  et  al .,  1998].  In  this  current  study,  the  CRRES  EFI 
data  are  used  to  determine  the  electric  power  spectral 
density  as  a  function  of  L  and  Kp.  These  spectra,  together 
with  various  assumptions  regarding  the  global  structure  of 
fluctuations,  are  used  to  estimate  the  radial  diffusion 
coefficient  Dll-  In  section  2,  the  CRRES  orbit  and  electric 
and  magnetic  field  instrumentation  are  briefly  reviewed. 
Section  3  describes  the  formalism  for  determining  Dll  in 
general,  and  outlines  the  approach  followed  here  necessi¬ 
tated  by  single  point  measurements.  The  determination  of 
the  power  spectral  densities  and  the  resulting  Dll  are 
provided  in  section  4.  Section  5  discusses  the  results,  and 
concluding  statements  are  made  in  section  6. 

2.  CRRES  Orbit  and  Instrument  Description 

[6]  CRRES  was  launched  on  July  25,  1990,  into  a 
geosynchronous  transfer  orbit  (GTO)  with  a  perigee  of 
350  km,  an  apogee  of  36000  km,  and  an  inclination  of 
18°.  With  an  orbital  period  of  10  hours,  it  made  a  near- 
equatorial  radial  pass  through  the  radiation  belts  about  four 
times  per  day.  The  CRRES  orbit’s  initial  apogee  at 
~0600LT  precessed  through  local  midnight  to  ~1000LT 
in  the  14  months  of  operation.  After  the  initial  spin-up 
period  following  launch,  the  satellite  spin  period  stabilized 
at  approximately  2  rpm  with  the  spacecraft  spin  axis  aligned 
to  within  15°  of  the  Earth-Sun  line. 

[7]  CRRES  carried  instrumentation  to  measure  both  elec¬ 
tric  and  magnetic  fields.  The  electric  field  instrument  (EFI) 
[Wygant  et  al 1992]  measured  the  electric  potential  be¬ 
tween  a  pair  of  spherical  sensors  and  between  a  pair  of 
cylindrical  antennas,  each  orthogonal  to  the  CRRES  spin 
axis.  The  EFI  signal  was  sampled  at  a  rate  of  32  samples/s 
and  filtered  at  1 0  Hz,  providing  a  measurement  of  the  two- 
dimensional  electric  field  (in  the  spacecraft  spin  plane)  with 
a  sensitivity  of  better  than  0. 1  mV/m  and  a  dynamic  range 
of  1000  mV/m.  The  full  resolution  data  was  spin-fit  to 
remove  the  sinusoidal  variation  of  the  electric  field  mea¬ 
surement  resulting  from  the  spacecraft  spin.  This  procedure 
yielded  the  electric  field  components  Ey  and  Ez  in  a 
modified  GSE  coordinate  system  (defined  in  section  4.1). 

[8]  The  fluxgate  magnetometer  [Singer  et  al. ,  1992]  was 
comprised  of  three  orthogonal  sensors  located  at  the  end  of 
a  6.1  m  boom,  and  yielded  2-s  resolution  of  the  full 
magnetic  field  vector.  The  magnetometer  operated  in  low 
(high)  gain  mode  with  a  dynamic  range  of  ±45000  nT 
(±850  nT)  for  CRRES  below  (above)  L  ~  3.5  RE.  The 
magnetic  field  vector  data  permitted  the  estimation  of  the 
component  of  the  electric  field  vector  out  of  the  spin  plane 
from  the  assumption  E||  =  0. 

3.  Radial  Diffusion  Coefficient  Formalism 

[9]  A  large  number  of  small  electrostatic  or  electromag¬ 
netic  impulses  randomly  perturbing  an  azimuthally  drifting 
particle  will  lead  to  a  net  change  in  its  radial  position  that 
can  be  modeled  as  a  diffusive  process.  Within  the  radial 
diffusion  formalism,  the  detailed  nature  of  these  impulses  is 
ideally  prescribed  by  the  electrostatic  (Dll)  or  electromag¬ 
netic  (Dll)  diffusion  coefficients.  In  practice,  a  number  of 


simplifying  assumptions  must  be  made  regarding  the  spatial 
and  temporal  variations  of  the  fields  involved.  In  this  paper 
we  are  mainly  interested  in  the  contribution  of  electrostatic 
fluctuations  to  the  radial  diffusion  process.  For  simplicity, 
the  derivation  of  Dll  has  traditionally  assumed  a  symmetric 
dipole  background  magnetic  field  and  has  considered  only 
the  azimuthal  component  of  E  arising  from  convection 
electric  field  variations  [j Falthammar ,  1965;  Schulz  and 
Lanzerotti ,  1974].  However,  more  recently  Elkington  et  al. 
[2003]  have  investigated  other  forms  for  DLl  assuming  a 
compressed  dipole  field  and  toroidal  and  poloidal  ULF 
wave  modes  in  addition  to  a  background  convection  electric 
field.  They  have  found  that  radial  diffusion  and  energization 
rates  increase  with  increasing  radial  distortion  of  the  mag¬ 
netic  field  and  increasing  convection  electric  field.  These 
rates  can  be  at  least  as  significant  as  those  for  the  traditional 
symmetric  mode  of  diffusion,  and  should  be  considered 
when  evaluating  diffusion  rates  during  periods  when  the 
magnetosphere  experiences  significant  radial  distortion 
from  the  impact  of  solar  wind  shocks  or  the  interaction 
with  high  speed  solar  wind  streams.  In  this  study,  we  follow 
the  approach  outlined  by  [. Falthammar ,  1965]  based  upon 
linearized  theory  which  describes  the  net  effect  of  randomly 
varying  electric  potential  fields  in  a  symmetric  (dipole) 
field,  which  violate  the  third  invariant  while  conserving 
the  first  and  second  invariants. 

3.1.  Electromagnetic  Versus  Electrostatic  Fluctuations 

[10]  Falthammar  [1965]  begins  the  derivation  of  the 
radial  diffusion  coefficient  with  the  first-order  equation  for 
the  time-dependent  radial  position  (r): 

dr  _  (  £«>(r0,<M)  £Vo  db(rQ,<\>,t)\ 
dt  \  Bo  3S0  54>  )^,+a 

where  E^  is  the  azimuthal  component  of  the  electric  field, 
B0  is  the  zero-order  magnetic  field  at  radial  position  rQ,  £20  is 
the  azimuthal  drift  frequency,  and  b  is  the  z-component  of 
the  fluctuating  magnetic  field.  The  azimuthal  component  E0 
may  be  expressed  by  the  Fourier  series: 

N 

E^(r0^,t)  =  ^£^(ro,/)cos(mt>  +  y(ro,0)  (2) 

/i=i 

[11]  For  the  case  of  electromagnetic  fluctuations,  both 
terms  in  the  brackets  of  (1)  are  retained;  E  has  a  nonzero 
curl,  and  VxE  =  —  ^  is  used  to  rewrite  (1)  in  terms  of  b 
alone.  Falthammar  [1965]  assumes  a  model  storm  magnetic 
field  with  azimuthally  symmetric  and  asymmetric  parts,  and 
writes  the  magnetic  radial  diffusion  coefficient  D^l  in  terms 
of  the  power  spectra  of  the  asymmetric  component  of  the 
magnetic  disturbance. 

[12]  For  the  case  of  purely  electrostatic  fluctuations,  the 
second  term  in  the  brackets  of  (1)  is  time-independent  and 
thus  averages  out,  leaving  only  the  electrostatic  field  E*. 
Falthammar  [1965]  derives  an  expression  for  DlL  in  terms 
of  the  summation  over  the  power  spectra  of  the  Fourier 
coefficients  E^fo,  t)  in  (2). 

3.2.  Computation  of  DEL 

[13]  No  attempt  has  been  made  here  to  identify  and 
separate  those  components  of  the  measured  electric  field 
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fluctuations  with  zero  curl  (electrostatic)  and  nonzero  curl 
(inductive).  In  what  follows,  we  pursue  the  computation  of 
Dll  for  electrostatic  fields  as  prescribed  by  Falthammar 
[1965]  and  as  followed  by  Holzworth  and  Mozer  [1979]. 

[14]  Holzworth  and  Mozer  [1979]  rewrote  Falthammar's 
[1965]  series  expansion  for  E^  (2)  as: 

00 

£0(Z,,<M)  =  ^£0>(Z.,<)cosn4)  +  £®  (I,0sinn<t>  (3) 

n=l 

where  cj)  =  is  the  azimuthal  position  of  the  particle 
drifting  with  frequency  f.  As  pointed  out  by  Holzworth 
and  Mozer,  the  determination  of  coefficients  E*AB(t)  up  to 
the  n  =  mth  term  requires  that  E(t)(L,  <(>,  t)  be  measured  at  a 
minimum  of  2m  azimuthal  positions  for  several  hours 
continuously.  At  each  time  t,  the  2m  values  of  E0(L,  (\> ,  t) 
are  Fourier  analyzed  with  respect  to  c|>  to  determine 
EonB(L,  t),  which  are  then  Fourier  analyzed  with  respect  to 
t  to  determine  the  power  spectral  density  P„,B.  The  total 
power  P  is  determined  by  summing  over  all  contributing 
modes: 

P(f,L,Kp)  =  '£)Pl(nf,L,Kp)  +PBM,L,Kp)  (4) 

n=  1 

The  radial  diffusion  coefficient  is  thus  given  as: 

=  - -~-P{f,L,Kp),  (5) 

where  Beq  is  the  equatorial  magnetic  field. 

[15]  Given  the  fact  of  single  spacecraft  measurements 

available  here,  there  are  a  number  of  significant  approx¬ 
imations  made  in  computing  Dll-  First,  we  begin  by 

assuming  that  only  the  n  =  1  coefficients  E£jB(L,  t)  are 

nonzero.  Thus  we  have 

E «  (£,  <t>.  0  =  El  (£,  0  cos  4>  +  £* (L,  l)  sin  <t>  (6) 

where  the  coefficients  E0  B(L,  t)  =  E0jB(L,  t)  are  in 
principle  determined  for  each  time  t  by  Fourier  analyzing 
(in  azimuth)  measurements  made  simultaneously  at  multiple 
(>2)  longitudes  c().  Although  these  measurements  do  not 
exist,  CRRES  measurements  do  yield  the  azimuthal  electric 
field  component: 

E0  =  -Ex  sin  <t>  -I-  Ey  cos  <(>  (7) 

for  each  position  (L,  4>)  and  time  t  along  the  CRRES 

trajectory  (with  <j>  =  0  at  local  noon).  If  we  equate  (6)  and  (7) 
and  make  the  assumption  that  the  measured  Ex  and  Ey 
are  independent  of  4>,  then  we  can  set  E^(L,  t)  =  Ey  and 
Eb(L,  t)  =  — Ex.  We  make  the  further  assumption  that 
for  a  given  FFT  interval,  the  satellite  remains  at  a  fixed  L 
(a  good  approximation  at  high  L,  but  less  so  at  low  L).  We 
thus  determine  the  total  power  spectral  density  P  in  (5) 
from 

P  =  Pf+PB  =Py+P„  (8) 


where  Px  and  Py  are  computed  from  the  Fourier  transform 
of  Ex  and  Ey,  respectively. 

4.  Power  Spectral  Density 

[16]  In  this  section  we  discuss  the  details  for  computing 
the  x  and  y  components  of  the  power  spectral  density  (solar 
magnetic  coordinates),  Px  and  Py,  of  the  near-equatorial 
electric  field  fluctuations  (0.2-15.9  mHz)  as  a  function  of  L 
(L  =  3  to  7)  and  magnetic  activity  (Kp  =  0  to  7). 

4.1.  Data  Preparation  and  Coverage 

[17]  We  begin  with  the  30  s  spin-fit  electric  field  vector 
components  in  the  spin  plane  (Ey,  Ez)  represented  in  a 
vehicle-based  nonrotating  coordinate  system  (modified 
geocentric  solar  ecliptic,  MGSE).  The  MGSE  coordinate 
system  is  defined  by  an  x-axis  that  points  in  the  positive  spin 
axis  direction  (9°  from  the  Sun),  a  y-axis  that  intersects  the 
ecliptic  and  spin  planes  and  points  toward  dusk,  and  a  z-axis 
that  completes  a  right-hand  rectangular  coordinate  system. 
The  measured  E  components  are  provided  in  the  moving 
spacecraft  reference  frame  and  must  be  transformed  to  a 
stationary  reference  frame  by  subtracting  the  (v  x  B)  field, 
where  v  is  the  spacecraft  velocity  and  B  is  the  locally 
measured  magnetic  field.  The  assumption  of  E  B  =  0  is 
used  to  determine  EX(MGSE)  from  the  measured  Ey(MGSE), 
EZ(MGSE),  and  B  vector  field.  The  vector  components 
of  E  in  the  MGSE  system  are  then  transformed  to  the 
solar  magnetic  (SM)  coordinate  system  defined  by  a  z-axis 
coincident  with  the  magnetic  dipole  axis,  an  x-axis 
perpendicular  to  the  z-axis  and  in  the  plane  containing 
the  z-axis  and  the  Sun-Earth  line,  and  a  y-axis  that 
completes  the  right-hand  rule.  The  spectral  analysis  is 
performed  on  the  time  series  for  EX(SM)  and  Ey(SM). 

[is]  It  was  necessary  to  first  identify  continuous  intervals 
that  met  stringent  “quality  control”  criteria,  and  this 
dramatically  reduced  the  volume  of  data  available  for 
analysis.  There  was  an  electrical  short  in  the  spherical 
probe  system  that  resulted  in  an  electric  field  measurement 
that  was  dependent  upon  the  plasma  density  and/or  space¬ 
craft  potential.  Also  affecting  the  spherical  probes  were 
single  event  upsets,  primarily  during  1990.  Because  of 
these  two  factors,  it  was  decided  to  use  only  the  cylindrical 
probes. 

[19]  The  cylindrical  probes  were  susceptible  to  techni¬ 
cal  difficulties  as  well.  During  most  of  1990  the  cylin¬ 
drical  probes  were  improperly  biased.  Additionally,  during 
1990  the  CRRES  apogee  was  in  the  postmidnight  sector 
where  spacecraft  charging  during  energetic  electron  injec¬ 
tions  often  degraded  EFI  probe  measurements.  This  study 
is  therefore  restricted  to  cylindrical  probe  data  from 
January  through  October  1991  when  probe  biasing  was 
optimal.  This  represents  an  immediate  reduction  of  the 
data  set  by  approximately  one  third,  eliminating  measure¬ 
ments  from  early  morning  local  times.  The  local  time 
coverage  runs  from  ~0300LT,  through  midnight,  and  up 
to  1200LT,  thus  precluding  the  prenoon  sector  (0600LT  to 
1200LT). 

[20]  There  were  three  other  factors  independent  of  the  EFI 
performance  that  further  reduced  the  data  set.  One  was  the 
routine  spacecraft  maneuvers  that  rendered  the  EFI  signal 
useless  and  required  the  elimination  of  the  10-hour  interval 


3  of  15 


A02214 


BRAUTIGAM  ET  AL.:  CRRES  ELECTRIC  FIELD  POWER  SPECTRA 


A02214 


following  any  attitude  adjustment.  A  second  was  the  uncer¬ 
tainty  of  spacecraft  attitude  and  velocity  that  introduced 
error  via  the  transformation  of  the  measured  E  to  a  frame 
corotating  with  the  Earth  [Wygant  et  al ,  1998].  This  error 
maximizes  at  low  altitude  where  the  spacecraft  speed  is 
largest  and  therefore  this  study’s  radial  coverage  is  con¬ 
strained  to  L  >  2.75.  A  third  factor  was  that  EFI  measure¬ 
ments  were  limited  to  the  two  E  components  in  the 
spacecraft  spin  plane  (Ey,  Ez,  mGSE)  and  thus  necessitated 
approximating  the  third  component  (Ex,  mGSE)  along  the 
spin  axis.  Using  the  assumption  E  •  B  =  0  leads  to  Ex  = 
— (ByEy  +  BZEZ)/BX  which  introduces  potentially  large 
uncertainties  in  Ex  for  small  Bx.  The  requirement  on  Bx 
for  reliable  estimates  of  Ex  is  given  by  Bx  >  0.5  (By  +  Bz)1/2; 
all  data  points  failing  this  criterion  were  rejected.  This  led  to 
more  data  being  rejected  throughout  ~  1 800LT  to  2000LT, 
over  a  broad  range  in  L,  than  over  the  ranges  1300LT  to 
1700LT  and  2000LT  to  2400LT. 

[21]  The  total  number  of  data  points  in  the  complete 
CRRES/EFI  database  is  1.075  million.  About  32%  of  these 
were  rejected  because  they  occurred  in  1990  when  the 
biasing  problem  was  most  severe,  30%  were  rejected 
because  of  the  criterion  on  Bx,  17%  were  rejected  because 
they  were  positioned  below  L  =  2.75,  and  1.3%  were 
rejected  because  of  the  large  perturbations  caused  by 
spacecraft  attitude  adjustments.  The  final  number  of  remain¬ 
ing  points  is  ~0.2 1 6  million,  which  represents  ~20%  of  the 
entire  database.  Unfortunately,  not  all  of  these  “good”  data 
were  used  because  intervals  with  a  minimum  number  of 
contiguous  points  were  required  to  yield  a  valid  power 
spectrum. 

4.2.  FFT  Analysis 

[22]  The  power  spectral  densities,  Pf  =  Py  and  Pf  = 
Px  (as  in  equation  (8))  are  determined  from  the  Fast 
Fourier  Transform  (FFT)  of  e£(L,  t)  =  Ey  and  EB(L,  t)  = 
— Ex,  respectively,  where  the  time  series  intervals  consist  of 
N  (=2n)  data  points  uniformly  spaced  in  time.  As  is 
standard  practice,  any  DC  offset  is  first  subtracted  and  a 
window  function  (wk)  then  applied  to  the  time  series 
before  the  FFT  is  evaluated  [Press  et  al .,  1992].  We  found 
that  various  tapered  windows  all  gave  essentially  the  same 
result,  and  have  chosen  the  Welch  window  for  the  final 
analysis. 

[23]  Unfortunately,  the  literature  on  power  spectra  can 
be  confusing  because  there  is  no  single,  universally  used 
convention  defining  the  relation  between  an  FFT  and  its 
power.  Falthammar  [1965]  uses  the  one-sided  power 
spectrum  convention  of  Rice  [1954],  which  we  adhere  to 
as  well.  One  may  write  the  power  spectral  density  in  terms 
of  the  discrete  Fourier  transform  (including  the  window 
function)  as 


KB(f)  = 


2 

WCsf 


N- 1  A 

£  e* p(-*2ir/fcA/)  , 

*=o 


(9) 


where  W  is  a  normalization  constant,  W  =  N 

k=0 

determined  from  the  window  function.  Defining  At  as  the 
sampling  time,  the  discrete  time  t  =  kAt,  the  duration  of  the 
FFT  interval  is  T  =  NAt,  and  the  discrete  frequency  step  is 
Af  =  1/T.  It  should  be  noted  that  since  the  e£,b  are  real, 


Table  1.  Drift  Frequencies  and  Energies  for  Fixed  L  and  First 
Invariant 


L 

=  100  MeV/G 

H  =  5000  MeV/G 

H  =  68000  MeV/G 

3.0 

6.5 

0.7  MeV,  0.55  mHz 
0.1  MeV,  0.22  mHz 

7.3  MeV,  4.3  mHz 
2.0  MeV,  2.9  mHz 

28.1  MeV,  15.9  mHz 
8.5  MeV,  10.8  mHz 

Pf ,B(f)  =  Pf,B(-f),  and  (9)  is  evaluated  for  positive  f  only; 
the  factor  of  2  in  (9)  accounts  for  equal  amounts  of  power  at 
both  positive  and  negative  f. 

[24]  There  is  a  trade-off  between  spatial  and  frequency 
resolution  that  requires  some  consideration.  The  set  of 
positive  frequencies  is  defined  by  {1/T,  2/T,  ...,  N/(2T)}. 
The  Nyquist  frequency  fN  =  N/(2T)  =  1/(2 At)  is  therefore 
fixed  at  16.95  mHz  given  that  the  average  sampling  time 
At  =  29.5s,  whereas  the  frequency  resolution  Af  =  l/(NAt) 
depends  upon  the  chosen  interval  size  (N).  For  N  =  32 
(T  =  944s),  Af  =  1.06  mHz;  for  N  =  64  (T  =  1888s),  Af  = 
0.53  mHz;  and  for  N  =  128  (T  =  3776s),  Af  =  0.265  mHz. 
Of  particular  interest  is  the  frequency  range  of  electric 
field  fluctuations  that  resonate  with  the  azimuthally  drift¬ 
ing  radiation  belt  particles.  Table  1  specifies  (over  a  range 
of  L,  energy,  and  corresponding  drift  frequency)  the 
radiation  belt  electron  population  with  which  this  study’s 
E  fluctuations  may  resonate.  These  electron  energies  range 
from  0.1  MeV  at  L  =  6.5  (0.22  mHz)  corresponding  to 
the  relatively  common  substorm  injected  population,  to 
28.1  MeV  at  L  =  3.0  (15.9  mHz)  corresponding  to  extremely 
rare  acceleration  events  [Li  et  al .,  1993].  The  main  concern 
here  is  to  specify  the  lower  energy  limit  where  the  drift 
frequency  is  limited  by  the  power  spectra  frequency 
resolution.  Although  the  lower  frequency  limit  would 
argue  for  the  N  =  128  case,  the  spatial  resolution  would 
be  greatly  compromised.  The  time  for  CRRES  to  traverse 
a  fixed  width  in  L  varies  with  L  such  that  within  the  time 
interval  for  N  =  128  (3776s),  CRRES  may  pass  from  L  = 
2.75  to  5.25,  or  from  L  =  5.25  to  6.75.  Our  targeted  spatial 
resolution  is  ~0.5  RE,  which  is  more  in  line  with  N  =  32 
(Af  =  1 .06  mHz).  We  ultimately  chose  to  rely  primarily  on 
the  N  =  128  case,  but  not  until  we  performed  the  complete 
analysis  for  all  3  cases  (N  =  32,  64,  128)  and  discovered 
that  the  L  dependence  on  the  power  spectral  density  was 
relatively  weak  which  meant  that  the  spatial  resolution  was 
not  a  critical  factor. 

[25]  Besides  the  spatial  resolution,  the  interval  size  also 
determines  the  total  number  of  FFT  time  series  available 
for  the  study.  After  imposing  the  EFI  data  “quality 
control”  criteria  discussed  earlier,  a  total  of  (1190,  1006, 
743)  valid  segments  of  continuous  data  were  found  for  the 
3  values  of  N  of  interest  (32,  64,  128),  with  the  total 
number  of  points  per  segment  ranging  from  N  (the 
minimum  required)  to  824.  It  is  a  standard  procedure  to 
perform  FFTs  on  overlapping  intervals.  The  spectral 
variance  is  minimized  if  the  intervals  are  offset  by  N/2 
[Press  et  al ,  1992],  but  we  use  intervals  offset  by  N/4  to 
greatly  decrease  the  number  of  “excess”  data  points  not 
utilized.  This  yields  a  total  of  19274  16-point  spectra, 
7988  32-point  spectra,  and  2751  64-point  spectra  for 
analysis. 

[26]  An  average  L  and  Kp  are  determined  for  each  FFT 
interval,  and  are  used  to  sort  the  power  spectra  accord¬ 
ingly.  The  distribution  of  spectra  over  the  range  of  L  bins 
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Table  2.  Number  of  Power  Spectra  in  L  and  Kp  Binsa 


L  * 

o 

ii 

CL. 

* 

Kp  =  1 

Kp  =  2 

Kp  =  3 

II 

c. 

II 

a. 

Kp  =  6 

Kp  >  7 

3.0 

11,  0,  0 

73,  0,  0 

86,  0,  0 

123,  0,  0 

75,  0,  0 

33,  0,  0 

18,  0,0 

21,0,0 

3.5 

27,  6,  0 

191,  70,  1 

213,  76,2 

290,  113,  0 

191,  70,  0 

77,  29,  0 

54,  26,0 

53,  22,  0 

4.0 

26,  11,  3 

196,  84,  27 

250,  115,  38 

326,  139,  56 

207,  98,  39 

94,  30,  11 

69,  33,  17 

66,  26,  10 

4.5 

26,  14,  5 

199,  88,  38 

305,  132,  46 

363,  160,  53 

250,  114,  43 

116,  51,  16 

82,  35,  11 

78,  34,  10 

5.0 

34,  10,  2 

261,  110,  40 

369,  160,  59 

418,  182,  71 

317,  131,48 

152,  62,  18 

82,  35,  13 

81,  32,  9 

5.5 

49,  18,  4 

347,  143,  57 

449,  187,  68 

525,  233,  88 

455,  191,  68 

189,  74,  29 

94,  37,  1 1 

65,  23,  10 

6.0 

50,  19,  7 

487,  209,  78 

665,  297,  116 

674,  289,  115 

619,  268,  95 

238,  103,  34 

90,  40,  17 

108,  42,  8 

6.5 

68,  30,  1 1 

852,  364,  117 

1165,  514,  205 

1347,  563,  203 

971,  409,  149 

337,  128,  36 

154,  60,  24 

136,  57,  22 

7.0 

44,  15,  2 

592,  262,  97 

801,  351,  131 

860,  363,  130 

441,  198,  70 

227,  89,  26 

103,  49,  17 

169,  65,  20 

“Numbers  given  within  Kp  columns  are  for  N  =  32,  64,  and  128  case,  respectively.  L  bins  are  0.5  Re  wide,  centered  on  table  values.  Kp  =  0  bin  includes 
values  (0,  0.3);  Kp  =  1  bin  includes  values  (0.6,  1.0,  1.3),  and  so  on  for  higher  Kp. 


(0.5  RE  wide)  and  Kp  bins  (integral  values)  are  summa¬ 
rized  in  Table  2.  Intervals  with  an  average  Kp  >  7  have 
been  combined  into  a  single  Kp  bin.  The  fact  that  the  L 
bins  are  0.5  RE  wide  does  not  mean  that  the  entire  FFT 


interval  fits  within  that  bin,  but  only  that  its  midpoint 
does;  for  the  case  of  N  =  128  the  FFT  interval  may 
actually  span  between  1.5  and  2.5  RE.  Spectra  for  L  <  2.75 
were  eliminated  because  of  the  error  introduced  by  attitude 


Doys=  56.38400-  56.47171 


Interval:  0-127 

Day:  56.38400  -  56.42768 
L  =  6.13 
Kp  =  3.6 

Fit:  P=  0.002*^-1.243) 
cc=  0.957 


0.1  1.0  10.0 
f  [mHz] 

Figure  1.  Top  two  panels  show  sample  256-point  time  series  segment  of  Ex  and  Ey  (mV/m)  from  day 
56,  1991.  FFTs  are  performed  on  the  128-point  interval  (points  0  to  127),  and  the  resulting  spectra  (Px,  Py, 
and  P  =  Px  +  Py)  are  shown  in  the  bottom  panel.  The  endpoints  (in  decimal  day),  the  average  L,  and  the 
average  Kp  are  shown  for  this  interval.  The  power  law  fit  to  P  is  listed,  with  a  correlation  coefficient 
(0.957)  indicating  an  excellent  fit. 
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Days=233. 48730- 233. 57428 


Interval:  64  -  191 
Day:  233.50913  -  233.55245 
L  =  5.98 
Kp  =  3.6 

Fit:  P=  0.238*fA(- 1 .103) 
cc=  0.039 


0.1  1.0  10.0 
f  [mHz] 

Figure  2.  Same  format  as  Figure  1,  but  for  a  segment  on  day  233,  1991.  FFTs  are  performed  on  the 
128-point  interval  (points  64  to  191).  The  power  law  fit  to  P  is  listed,  with  a  correlation  coefficient 
(0.039)  indicating  a  poor  fit  due  to  the  strong  peak  at  2.7  mHz. 


uncertainties;  those  for  L  >  7.25  were  eliminated  because 
of  poor  statistics  in  that  L-bin. 

4.3.  Power  Spectral  Density  Model 

[27]  In  choosing  to  create  an  empirical  model  of  the 
electric  field  power  spectral  density,  and  hence,  the  radial 
diffusion  coefficient  Dll,  one  must  evaluate  how  best  to 
parameterize  the  power.  The  temporal  variability  may  in 
principle  be  parameterized  by  any  number  of  parameters. 
For  example,  Li  et  al.  [2001]  have  used  solar  wind  speed 
(Vsw)  fluctuations  to  derive  a  time-dependent  DLl,  and 
Mathie  and  Mann  [2001]  have  parameterized  ULF  power 
with  Vsw.  The  z-component  of  the  interplanetary  magnetic 
field  (IMF  Bz)  controls  the  magnitude  of  the  convection 
electric  field  and  may  thus  also  be  expected  to  determine  its 
fluctuation  power.  However,  throughout  the  CRRES  mis¬ 
sion  IMP-8  was  the  only  source  of  Vsw  and  IMF  Bz 
measurements  and  these  were  too  sparse  for  any  systematic 
studies.  The  global  magnetic  activity  index  Kp  has  been 


used  in  the  past  to  parameterize  electric  field  measurements 
[. Mozer ,  1971;  Rowland  and  Wygant ,  1998],  and  is  used  here 
as  well.  There  are  two  alternate  approaches  possible  to 
compute  the  model  radial  diffusion  coefficient  DlL(L,  Kp). 
The  first  is  to  calculate  the  power  spectra  and  the  respective 
Dll  for  each  FFT  time  interval,  and  then  sort  and  average 
the  individual  Dll’s  according  to  L  and  Kp.  The  second 
approach  is  to  first  create  an  average  power  spectra  model 
which  can  then  be  used  to  calculate  the  model  Dll-  We 
follow  the  latter  approach  because  the  model  power  spectra 
are  of  interest  in  their  own  right. 

[28]  Figures  1  and  2  show  representative  Ex  and  Ey  time 
series  and  their  related  power  spectra.  The  signal  in  Figure  1 
(day  56)  has  a  significant  DC  offset,  particularly  obvious  in 
the  first  half  of  the  256-point  segment  (sample  points  0- 
127)  for  the  Ex  component.  The  maximum  peak-to-peak 
amplitude  Ex  and  Ey  increases  from  less  than  0.5  mV/m 
during  the  first  half  of  segment,  to  1-1.5  mV/m  during  the 
second  half.  The  summed  power  spectrum  (Px  plus  Py)  from 
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Scaled  log(Py) 


Figure  3.  Histogram  of  the  logarithmic  power  (Py)  at 
2.65  mHz  for  the  model  bin  (L  =  6.0,  Kp  =  4).  The 
lognormal  distribution  defined  by  a  mean  equal  to  —1.76, 
and  standard  deviation  equal  to  0.85  is  superimposed. 


the  N  =  128  interval  (data  points  0-127)  is  well  fit  by  a 
power  law  in  frequency  with  an  index  of  —1.2  and  com¬ 
pares  well  with  the  average  power  law  index  of  -1.6  ±  0.3 
found  by  Mozer  [1971],  This  time  interval  has  an  average 
L  (=  6.13)  and  Kp  (=  3.6)  that  assigns  this  spectrum  to 
the  L  =  6.0,  Kp  =  4  model  bin.  The  signal  in  Figure  2 
(day  233),  assigned  to  the  same  model  bin  but  separated 
by  177  days,  shows  very  different  characteristics.  Its 
maximum  peak-to-peak  value  is  4-6  mV/m,  significantly 
larger  than  the  signal  in  Figure  1.  Furthermore,  there  is  a 
coherent  oscillation  evident  in  the  signal  that  leads  to  a 
spectral  peak  around  2.8  mHz,  thus  rendering  a  poor  fit  to 
a  power  law.  These  coherent  oscillations  (with  a  spectral 
peak  between  2  and  10  mHz)  are  common  enough  that  the 
model  average  spectra  exhibit  a  broad  peak  superimposed 
on  an  otherwise  power  law  form.  Rather  than  attempting 
to  derive  an  analytic  expression  for  this  spectral  shape, 
we  fit  the  model  power  at  each  frequency  as  a  function 
of  (L,  Kp). 

[29]  For  each  case  of  N  (32,  64,  128),  a  model  spectrum 
for  Px(f)  and  Py(f)  is  computed  for  each  (L,  Kp)  bin  by 
averaging  together  the  logarithm  of  the  power  spectra 
assigned  to  that  given  bin.  The  number  of  spectra  in  each 
(L,  Kp)  bin,  for  each  N,  is  summarized  in  Table  2.  The 
model  power  spectrum  for  (L  =  6.0,  Kp  =  4)  is  constructed 
from  95  individual  spectra  (for  N  =  128);  and  the  histogram 
for  log  Py(f  =  2.65  mHz)  is  shown  in  Figure  3.  The  data  are 
well  approximated  by  the  superimposed  normal  distribution 
defined  by  a  mean  of  —  1 .76  and  standard  deviation  (sigma) 
of  0.85.  The  model  spectrum  Py(L  =  6.0,  Kp  =  4)  for  N  = 
128  is  shown  in  Figure  4.  The  standard  deviation  from  the 
average  power  at  a  given  frequency  is  represented  by 
vertical  bars,  and  is  seen  to  increase  slightly  with  frequency. 
The  log  power  at  all  frequencies  for  all  model  spectra  is 


characterized  by  a  sigma  on  the  order  of  1 ,  indicating  that 
the  observed  power  is  often  (~16%  of  the  time)  up  to  a 
factor  of  10  times  higher  (or  lower)  than  that  represented  by 
the  model. 

[30]  As  discussed  in  section  4.2,  it  is  preferable  to  utilize 
the  N  =  128  case  because  its  corresponding  lower  frequency 
limit  provides  better  coverage  for  the  ~  1 00  keV  electrons  of 
interest  to  radiation  belt  dynamics.  However,  as  indicated 
in  Table  2,  the  N  =  128  case  yields  essentially  no  spectral 
data  in  the  2  lowest  L  bins  (L  =  3.0  and  3.5).  Likewise, 
the  N  =  64  case  yields  no  data  in  the  L  =  3.0  bin.  To 
rectify  this  situation,  the  N  =  64  (32)  case  is  used  to 
provide  spectral  data  for  the  L  =  3.5  (3.0)  bin,  and  the 
available  spectra  are  extrapolated  to  the  desired  lower 
frequency  limit  (0.20  mHz).  For  the  N  =  32  case,  the 
power  at  the  first  3  frequencies  (1.06,  2.12,  and  3.18  mHz) 
is  fit  to  a  power  law  (usually  a  good  approximation  below 
3.18  mHz)  and  extrapolated  to  0.53  and  0.20  mHz. 
Likewise,  for  the  N  =  64  case,  the  power  at  the  first  4 
frequencies  (0.53  to  2.12  mHz)  is  fit  to  a  power  law  and 
extrapolated  to  0.20  mHz.  The  spectra  for  the  N  =  128 
case  have  a  lower  frequency  of  0.26  mHz  and  to  establish 
a  consistent  set  of  frequencies,  the  power  between  0.26 
and  2.12  mHz  is  fit  to  a  power  law  and  extrapolated  to 
0.20  mHz. 

[31]  To  examine  the  validity  of  the  extrapolation  to  the 
lowest  frequencies  using  N  =  32  spectra  (at  L  =  3.0)  and 
N  =  64  spectra  (at  L  =  3.5),  the  extrapolated  values  were 
directly  compared  to  the  N  =  128  spectra  case  for  L  =  4. 
This  comparison  showed  that  the  extrapolation  to  0.2  and 
0.53  mHz  at  L  =  3.0  underestimates  the  spectral  density 
by  up  to  a  factor  of  ~14  and  ~8,  respectively,  and  that 
the  extrapolation  to  0.2  mHz  at  L  =  3.5  underestimates 
the  spectral  density  by  up  to  a  factor  of  ~6. 


0.1  1.0  10.0 
f  [mHz] 

Figure  4.  Model  spectrum  (Py)  for  L  =  6.0  and  Kp  =  4 
with  bars  representing  the  standard  deviation  from  model 
average. 
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Figure  5.  Model  spectra  P  =  Px  +  Py  for  Kp  =  1,  3,  and  6  at  L  =  3.0  to  6.5  in  steps  of  0.5  RE.  The 
symbols  represent  the  model  averages;  the  lines  represent  the  analytic  fit  to  model  averages,  P  = 
aLbexp(cKp),  where  each  frequency  is  independently  fit  as  a  function  of  L  and  Kp  (i.e.,  the  lines  are  not 
fits  of  the  form  power  versus  frequency). 


[32]  Considering  the  slow  variation  with  frequency 
exhibited  by  the  model  spectra,  it  is  felt  that  maintaining 
a  resolution  of  0.26  mHz  for  the  entire  frequency  range 
is  unwarranted.  Furthermore,  since  the  variance  in  power 
is  much  greater  at  the  Nyquist  frequency  than  below,  we 
restrict  the  study  to  below  the  Nyquist  frequency.  The 
frequencies  used  in  the  model  are  chosen  to  be  the  set 
in  common  (after  being  extrapolated)  with  the  three 
cases  of  N,  and  are  as  follows:  0.20,  0.53,  1.06,  2.12, 
...  15.9  mHz,  where  a  1.06  mHz  resolution  is  main¬ 
tained  above  1.06  mHz.  Examples  of  the  model  total 
power  spectra  P  =  Px  +  Py  for  various  (L,  Kp)  are 


illustrated  in  Figure  5.  As  noted  earlier,  the  spectra  deviate 
from  a  power  law  above  ~2-3  mHz  as  a  result  of  the 
spectral  peak  from  coherent  oscillations  in  the  electric 
field. 

4.4.  Power  Spectral  Density  Fits 

[33]  The  dependence  of  the  total  power  P  on  (L,  Kp) 
evident  in  Figure  5  is  more  closely  examined  in  Figure  6  at 
three  selected  frequencies.  Whereas  the  power  at  any  given 
frequency  may  vary  by  a  factor  of  2  to  10  over  the  full  range 
of  L  (3.0-7.0)  for  fixed  Kp,  it  may  vary  by  more  than  two 
orders  of  magnitude  over  the  full  range  of  Kp  (0  to  7)  for 
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Figure  6.  Plots  of  P  =  Px  +  Py  (at  0.2,  2.12,  and  8.47  mHz)  versus  (left)  L  for  (top  to  bottom)  Kp  =  1,  2, 
4,  and  6;  and  (right)  P  versus  Kp  for  (top  to  bottom)  L  =  3.5,  4.5,  5.5,  and  6.5.  The  symbols  represent  the 
model  averages;  the  lines  represent  the  analytic  fit  to  model  averages,  P  =  aLbexp(cKp),  where  each 


frequency  is  independently  fit  (as  in  Figure  5). 


fixed  L.  There  is  least  change  in  P  at  0.2  mHz  when  viewed 
as  a  function  of  L. 

[34]  In  fitting  P  as  a  function  of  (L,  Kp),  the  range  of 
(L,  Kp)  is  limited  to  L  =  3.5  to  6.5  and  Kp  =  1  to  6  where 
optimal  data  exist.  The  bins  with  Kp  =  0  are  avoided 
because  of  low  statistics,  and  the  Kp  >  7  bins  are  avoided 
because  they  represent  a  composite  bin  of  Kp  =  7,  8,  9. 
The  L  =  3.0  bins  are  avoided  because  of  the  uncertainty 
associated  with  the  (v  x  B)  correction  (greatest  below 
L  =  2.75)  compounded  by  relatively  low  statistics,  and 
L  =  7.0  bins  are  avoided  because  of  the  uncertainty 


related  to  off-equatorial  effects.  The  total  power  at  each 
frequency  is  fit  to  the  function  P(L,  Kp)  =  a  Lb  exp(cKp) 
using  a  singular  value  decomposition  algorithm  [Press  et 
al.,  1992].  We  were  first  motivated  to  use  the  exponential 
dependence  on  Kp  by  Mozer  [1971],  who  used  P(f,  Kp)  = 
a  exp(c  Kp)  f n,  and  we  found  that  this  did  an  adequate 
job  of  fitting  the  data  (right-hand  panels  of  Figure  6). 
Unlike  Mozer  [1971],  we  did  observe  a  slight  dependence 
on  L,  and  we  found  that  the  Lb  dependence  was  a  simple 
form  that  gave  a  reasonably  good  fit  to  the  data  (left-hand 
panels  of  Figure  6).  This  set  of  parameters  {a,  b,  c}  for 
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Table  3.  Fit  Parameters  for  Fixed  Frequency 


Fit  Parameters  of  P  =  a  Lb  exp(cKp) 


Frequency,  mHz 

a  [(mV/m)2/mHzl 

b 

c 

Corn  Coeff. 

0.20 

1.09  • 

KT1 

-1.31 

0.39 

0.81 

0.53 

1.56  • 

1(T3 

0.42 

0.53 

0.89 

1.06 

1.08  • 

IO"4 

1.53 

0.60 

0.79 

2.12 

3.24  • 

1(T6 

3.26 

0.71 

0.90 

3.18 

1.14  • 

10~6 

3.77 

0.76 

0.80 

4.24 

1.35  • 

10"6 

3.38 

0.86 

0.74 

5.30 

2.10  • 

10-6 

2.92 

0.86 

0.79 

6.36 

5.21  • 

10-6 

2.23 

0.85 

0.67 

7.42 

6.71  • 

1(T6 

1.85 

0.89 

0.81 

8.47 

8.45  • 

I0“6 

1.73 

0.83 

0.76 

9.53 

8.97  • 

10“6 

1.54 

0.87 

0.84 

10.6 

1.06  • 

1 0‘5 

1.38 

0.84 

0.76 

11.7 

5.18  • 

10-6 

1.56 

0.91 

0.86 

12.7 

4.39  • 

10"6 

1.57 

0.93 

0.78 

13.8 

6.12  • 

10“6 

1.38 

0.88 

0.87 

14.8 

3.91  • 

IO"6 

1.52 

0.94 

0.81 

15.9 

2.64  • 

10“6 

1.67 

0.94 

0.85 

each  model  frequency  is  provided  in  Table  3  and  is  plotted 
in  Figure  7.  There  is  a  steep  gradient  in  the  parameters 
(increase  in  b  and  c,  decrease  in  a)  up  to  3-4  mHz,  and 
above  ~6  mHz  they  maintain  a  fairly  constant  value.  This 
turning  point  at  3-4  mHz  coincides  with  the  point  at 
which  the  model  spectra  in  Figure  5  diverge  from  the 
power  law  behavior  observed  at  lower  frequencies  as  a 
result  of  the  type  of  peaked  spectra  illustrated  in  Figure  2. 
A  comparison  between  the  model  and  fitted  values  is 
shown  in  Figures  5  and  6.  In  Figure  6  the  comparison 
extends  throughout  the  full  range  of  L  and  Kp.  Within  the 
range  of  and  L  and  Kp  used  in  the  fitting  procedure  the 
agreement  is  generally  within  a  factor  of  2. 

[35]  Figure  8  illustrates  the  degree  to  which  this  dynamic 
(Kp-driven)  model  reproduces  the  observed  power  for  a 
given  frequency  (2.12  mHz)  and  L  (6.0).  To  quantify  the 
comparison  we  compute  the  following  quantity  as  a  function 
of  the  fitted  model  power  (Pflt)  and  the  power  computed  at 
each  of  the  ith  FFT  intervals  (Pdata)  for  a  given  L  and 
frequency,  with  the  sum  over  the  entire  period  of  1 99 1  under 
study: 


(1  fl/rts  2\ ^ 

—  -  log^,)  J  (io) 


the  electric  fields  were  not  mapped  to  the  equator  (see 
section  5),  the  computed  P  are  not  strictly  equatorial,  but  are 
up  to  20°  off  the  magnetic  equator.  Therefore,  for  consis¬ 
tency  the  Dll  is  determined  not  with  Beq  but  with  an 
average  CRRES  measured  B,  BCrres>  determined  as  a 
function  of  L  independent  of  Kp.  The  fitted  expression  for 
Bcrres  is  also  dipolar,  with  BCrres  =  4.56  •  104  L-3  [nT]. 

[37]  Radial  diffusion  coefficients  Dll(£  L,  Kp)  for  Kp  =  1 , 
3,  and  6  are  plotted  versus  L  in  Figure  9  for  first  adiabatic 
invariants  p  =  500  MeV/G  (top)  and  p  =  5000  MeV/G 
(bottom).  For  a  given  p,  each  value  of  L  is  associated 
with  an  equatorially  mirroring  particle  of  specified  energy 
and  corresponding  drift  frequency  given  by  [Schulz  and 
Lanzerotti,  1974]: 


(ii) 


where  y  is  the  relativistic  factor,  c  is  the  speed  of  light, 
Re  is  the  Earth  radius,  m0  is  the  electron  rest  mass,  q  is 


For  the  case  plotted  in  Figure  8,  6modej  =  0.64.  Consistent 
with  the  standard  deviations  discussed  in  the  preceding 
section,  a  moderate  amount  of  discrepancy  (up  to  a  factor 
of-10  too  low  or  high)  is  seen  between  Pdata  and  Pflt.  The 
values  of  6mode|  for  all  (f,  L)  are  given  in  Table  4.  In 
constructing  this  table’s  entries  for  L  =  3.0  and  L  =  3.5, 
where  32-point  and  64-point  intervals  are  used  respectively, 
the  spectra  are  not  extrapolated  to  the  lower  frequencies 
associated  with  the  128-point  intervals;  hence,  the 
corresponding  “no  data”  table  entries.  This  table  shows 
that  6modei  increases  with  frequency;  a  trend  that  is  consis¬ 
tent  with  the  larger  standard  deviations  observed  at  higher 
frequencies  as  noted  in  Figure  4. 

4.5.  Radial  Diffusion  Coefficients 

[36]  With  P  now  determined,  we  can  calculate  the 
radial  diffusion  coefficient  Dll  using  equation  (5).  Since 


Figure  7.  Fit  parameters  (top)  a,  (middle)  b,  and  (bottom) 
c  used  in  the  fit  P  =  aLbexp(cKp),  plotted  versus  frequency. 
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L  =  6.0  2.12  mHz 


Figure  8.  Comparison  of  observed  power  versus  fitted 
model  power  P(f  =2.12  mHz,  L  =  6.0)  for  the  full  period  of 
1991  under  study.  The  power  observed  for  individual  FFT 
intervals  (Pdata)  is  plotted  as  points,  and  the  Kp-driven 
model  power  (Pf,t)  is  plotted  as  a  solid  line.  The  value  of 
Kp  used  to  drive  the  model  is  that  averaged  over 
individual  FFT  intervals. 

the  electron  charge,  and  B0  is  the  equatorial  magnetic 
field  at  1  Re.  The  energies  (frequencies)  corresponding  to 
500  MeV/G  range  from  0.4  MeV  (0.8  mHz)  at  L  =  6.5  to 
2.0  MeV  (1.3  mHz)  at  L  =  3.0.  Likewise,  the  energies 
(frequencies)  corresponding  to  5000  MeV/G  range  from 
2  MeV  (2.9  mHz)  at  L  =  6.5,  to  7.3  MeV  (4.3  mHz)  at 
L  =  3.0.  It  should  be  recalled  that  the  values  of  DEl 
computed  for  particle  drift  frequencies  f  <  1  mHz  at  L  =  3.0 
(3.5)  rely  on  extrapolated  power  spectral  densities  from 
the  N  =  32  (64)  case;  otherwise,  the  N  =  128  case  is 
used.  For  the  lowest  frequencies,  the  power  (and  Dll) 


may  be  underestimated  by  as  much  as  an  order  of 
magnitude. 

5.  Discussion 

[38]  The  procedure  for  mapping  a  purely  potential  con¬ 
vection  electric  field  along  an  equipotential  magnetic  field 
line  between  the  ionosphere  (i.e.,  higher  latitudes)  and  the 
equator  has  been  worked  out  by  Mozer  [1970].  Earlier  work 
on  electric  power  spectra  [Mozer,  1971;  Holzworth  and 
Mozer ,  1979]  involved  balloon  measurements  that  required 
implementing  this  mapping  procedure  to  infer  the  power  at 
the  equator.  Since  CRRES  was  near  the  magnetic  equatorial 
plane  (between  0°  and  20°  magnetic  latitude),  and  because 
of  the  uncertainties  regarding  the  assumption  of  equipoten¬ 
tial  field  lines,  no  mapping  was  performed  for  this  study. 
Unlike  the  L  dependence  of  P  observed  in  this  study,  Mozer 
[1971]  inferred  from  his  ionospheric  measurements  that  P, 
(singly  measured  ith  component)  at  the  equator  was  inde¬ 
pendent  of  L  (for  L  =  3  to  8).  He  also  found  that  by 
comparing  power  at  local  noon  and  midnight  that  P,  is 
independent  of  local  time,  something  that  we  have  not 
examined  in  this  study.  As  noted  in  section  4.1,  there  is  a 
significant  deficit  of  LT  coverage  for  various  reasons.  In 
particular,  there  is  essentially  no  coverage  on  the  dawn  side 
from  0300LT  to  1 200LT,  and  minimal  coverage  on  the  dusk 
side  from  1 800LT  to  2000LT.  It  has  been  found  from  space- 
based  measurements  that  the  predominance  of  ULF  wave 
power  is  found  in  the  0600LT  to  1200LT  sector  [Anderson 
et  al. ,  1990],  while  ground  measurements  suggest  that 
during  moderately  active  magnetic  conditions  (-100  nT  < 
Dst  <  -50  nT)  ULF  wave  power  is  at  relatively  low  levels 
for  the  6  hours  centered  on  noon  [O'Brien  et  al .,  2003,  Plate 
2].  The  power  and  diffusion  coefficients  reported  in  this 
paper  may  be  lower  than  those  if  all  local  times  had  been 
sampled,  but  it  is  difficult  to  quantify  the  degree  to  which 
the  results  may  be  affected. 

[39]  In  section  4.2,  the  trade-off  between  higher  spatial 
resolution  (gained  by  shorter  FFT  intervals)  versus  higher 
frequency  resolution  (gained  by  longer  FFT  intervals)  was 
discussed.  In  Figure  10  we  plot  the  value  of  Py(f  =  2. 1  mHz) 


Table  4.  Values  of  6modci  (Equation  (10)) 


f,  Hz 

L  -  3.0“ 

L  =  3.5b 

L  -  4.0C 

L  =  4.5 

0 

11 

L  =  5.5 

L  -  6.0 

L  *■  6.5 

L  =  7.0 

0.20 

no  data 

no  data 

0.51 

0.58 

0.57 

0.65 

0.61 

0.62 

0.70 

0.53 

no  data 

0.51 

0.58 

0.61 

0.57 

0.60 

0.56 

0.59 

0.75 

1.06 

0.64 

0.54 

0.55 

0.52 

0.56 

0.60 

0.58 

0.58 

0.67 

2.12 

0.60 

0.59 

0.56 

0.60 

0.62 

0.67 

0.64 

0.64 

0.70 

3.18 

0.59 

0.57 

0.57 

0.59 

0.63 

0.71 

0.72 

0.68 

0.70 

4.24 

0.61 

0.62 

0.61 

0.69 

0.76 

0.72 

0.68 

0.73 

0.81 

5.30 

0.62 

0.68 

0.75 

0.71 

0.71 

0.79 

0.81 

0.73 

0.81 

6.36 

0.67 

0.65 

0.72 

0.74 

0.74 

0.84 

0.80 

0.75 

0.84 

7.42 

0.68 

0.72 

0.75 

0.79 

0.81 

0.82 

0.88 

0.80 

0.82 

8.47 

0.67 

0.67 

0.75 

0.71 

0.82 

0.83 

0.85 

0.84 

0.91 

9.53 

0.74 

0.75 

0.72 

0.74 

0.84 

0.90 

0.84 

0.86 

0.89 

10.6 

0.75 

0.77 

0.76 

0.76 

0.89 

0.94 

0.85 

0.87 

0.92 

11.7 

0.77 

0.82 

0.80 

0.84 

0.86 

0.85 

0.87 

0.85 

0.93 

12.7 

0.84 

0.77 

0.76 

0.80 

0.86 

0.91 

0.91 

0.83 

0.94 

13.8 

0.86 

0.81 

0.81 

0.78 

0.85 

0.88 

0.86 

0.83 

0.93 

14.8 

0.89 

0.78 

0.82 

0.75 

0.87 

0.88 

0.88 

0.88 

0.95 

15.9 

0.88 

0.79 

0.84 

0.82 

0.89 

0.94 

0.87 

0.85 

0.98 

“L  =  3.0  statistics  performed  with  N  =  32  intervals. 
bL  =  3.5  statistics  performed  with  N  =  64  intervals. 
CL  >  4.0  statistics  performed  with  N  =  128  intervals. 
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Figure  9.  Plots  of  Dll  for  fixed  first  adiabatic  invariant 
(top)  \i  =  500  MeV/G  and  (bottom)  p,  =  5000  MeV/G  as  a 
function  of  L,  for  Kp  =  1,3,  and  6.  The  symbols  (keyed 
to  Kp)  with  the  solid  lines  drawn  through  them  represent  the 
Dll  computed  with  the  model  averages;  the  dashed  lines 
represent  the  Dll  computed  with  the  fitted  model  using  the 
parameters  in  Table  3  (i.e.,  the  dashed  lines  are  not  the  “best 
fit”  line  through  the  symbols). 


determined  by  a  model  based  on  N  =  128  versus  that  for 
N  =  32,  at  three  different  L  intervals.  Although  the  power 
for  N  =  128  is  consistently  greater  than  that  for  N  =  32 
(by  less  than  a  factor  of  2),  there  is  no  significant  bias  in 
L  that  would  affect  the  L  dependence.  Thus  it  appears 
that  little  was  sacrificed  in  opting  for  the  greater  frequency 
resolution. 

[40]  Mozer  [1971]  found  that  the  power  in  a  single  ith 
component  of  electric  field  was  given  by  Pj(f,  Kp)  = 
Poexp(0.4Kp)  •  rn  with  PG  =  1.28  •  10"2  (mV/m)2/mHz, 
n  =  1.6,  and  with  f  in  mHz.  As  discussed  in  section  4.3  and 
illustrated  in  Figures  2  and  5,  spectra  for  moderate  to  high 


Kp  exhibit  peaks  in  the  several  mHz  range  and  result  in 
poor  power  law  fits.  These  peaks  are  generally  much  less 
pronounced  in  the  Kp  =  1  spectra  which  also  show  little 
variation  in  L.  For  comparison  our  Kp  =  1  model  spectra 
(Py)  are  averaged  over  all  L  to  give  Py  =  1.52  •  10-3f-M; 
Mozer ’s  result  (evaluated  for  Kp  =  1)  gives  Pj  =  1.91  • 
10”3Pl  6.  At  0.2  mHz,  Mozer ’s  result  is  ~25  times  larger 
than  our  result,  with  this  factor  decreasing  to  ~5  at  16  mHz. 

[41]  The  conversion  of  the  power  spectra  to  a  radial 
diffusion  coefficient  Dll(£  L,  Kp)  is  much  more  problem¬ 
atic  than  the  computation  of  the  power  spectra,  and  is 
difficult  to  meaningfully  quantify  the  errors  involved.  There 
are  two  major  sources  of  uncertainty,  and  we  merely  raise 
the  issues  here  rather  than  attempt  to  resolve  them.  One 
uncertainty  arises  from  the  complete  ignorance  of  the 
azimuthal  dependence  of  the  fluctuations  as  discussed 
above.  Lanzerotti  and  Wolfe  [1980]  point  out  that  single 
point  measurements  of  magnetic  or  electric  field  fluctua¬ 
tions  will  lead  to  an  overestimate  of  radial  diffusion  rates 
because  a  certain  fraction  of  power  at  the  particle  drift 
frequency  will  reside  in  higher  frequency  spatial  compo¬ 
nents.  There  are  alternatives  to  the  azimuthal  dependence 
assumed  in  this  study,  and  it  is  not  clear  how  these  various 
models  would  affect  the  resulting  Dll-  The  Volland-Stem 
convection  electric  field  model  [ Volland ,  1978]  describes  a 
field  with  a  component  in  the  dawn-dusk  direction  only, 
which  CRRES  data  clearly  show  is  not  the  case.  Riley  and 
Wolf[  1992]  provide  an  alternative  field  model  in  terms  of  a 
finite  Fourier  expansion  in  azimuthal  angle.  Another  model 
for  defining  the  azimuthal  dependence  could  be  derived 
through  the  statistical  analysis  of  MHD  simulations.  The 
approach  followed  in  this  study  is  only  one  way  of  dealing 
with  the  severe  handicap  of  single  point  measurements. 
A  second  uncertainty  is  the  determination  of  the  relative 
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Figure  10.  Plot  of  Py  at  2.1  mHz  computed  from  an  FFT 
interval  with  N  =  128  (Py28)  versus  that  from  an  interval 
with  N  =  32  (Py2).  The  comparison  is  done  for  three  L 
intervals  (L  =  4.5,  5.5,  and  6.5)  as  designated  by  the  legend 
key.  The  solid  line  corresponds  to  a  ratio  of  1  and  is 
included  to  aid  the  eye  in  comparison. 
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Figure  11.  DLL  for  fixed  first  adiabatic  invariant  p  =  100,  500,  1000,  and  5000  MeV/G.  Plotted  are 
the  Dll  [CRRES]  from  this  study  (solid  lines)  as  compared  to  the  DLlE  [B-A]  (dashed  lines)  and 
Dll  [B-A]  (dotted  lines)  defined  by  Brautigam  and  Albert  [2000],  for  Kp  =  1  (square)  and  Kp  =  6 
(triangle).  The  Dll  [CRJIES]  are  drawn  as  best  fits  to  the  form  DLlE  =  D0Ln. 


contribution  from  electrostatic  and  electromagnetic  fluctua¬ 
tions.  Progress  towards  quantifying  these  uncertainties  must 
come  from  future  analysis  between  electric  and  magnetic 
field  data. 

[42]  Figure  1 1  compares  this  study’s  results  Dll  [CRRES] 
with  the  Kp-dependent  coefficients  used  previously  by 
Brautigam  and  Albert  [2000]  (hereinafter  referred  to  as 
B-A)  to  model  the  dynamics  of  an  isolated  magnetic 
storm.  At  p  =  100  MeV/G,  for  both  Kp  values  shown, 
Dll  [B-A]  is  comparable  to  or  exceeds  DlL  [B-A].  As  p 
increases,  Dm  [B-A]  decreases  and  therefore  its  magnitude 
relative  to  Dll  (which  is  independent  of  p)  decreases.  For 
all  values  of  p  shown  (for  Kp  =  1),  Dll  [CRRES]  exceeds 
Dll  [B-A].  However,  Dll  [B-A]  exhibits  a  stronger 
dependence  on  Kp  so  that  for  Kp  =  6,  DlL  [B-A]  exceeds 
d(l  [CRRES]  for  all  but  p  =  5000  MeV/G. 

[43]  In  Brautigam  and  Albert  [2000]  the  case  was  made 
using  DLl  [B-A]  =  Dll  [B-A]  +  DlL  [B-A]  that  the  modeled 
radial  diffusion  could  account  for  the  p  =  100  MeV/G 


dynamics  reasonably  well,  but  could  not  explain  the  dynam¬ 
ics  (formation  of  peak  in  phase  space  density  at  L  ~  4.5)  for 
p  =  1000  MeV/G.  As  argued  below,  using  Dll  [CRRES]  in 
place  of  D[l  [B-A]  is  not  likely  to  affect  the  general 
conclusions  of  that  study.  It  is  primarily  the  total  DLl  at 
the  higher  activity  levels  (Kp  ~  6)  that  affects  the  final  result, 
so  for  p  =  100  MeV/G  (Kp  =  6),  where  Dll  [CRRES] 
Dll  [B-A]  »  D&  [B-A],  the  total  DLL  would  be  ~l/2 
that  of  the  original  study,  which  would  not  affect  the 
qualitative  result  for  low  p.  For  p  =  1000  MeV/G  (Kp  = 
6)  where  Dll  [CRRES]  <  D[l  [B-A]  <  Dll  [B-A],  the 
total  DLl  would  be  comparable  to  that  of  the  original 
study  (being  dominated  by  DlL)  and  would  thus  be  not 
likely  to  affect  the  results  for  high  p. 

[44]  In  addition  to  past  estimates  of  Dll  and  Dll  from 
field  measurements,  there  have  been  estimates  of  DLl  from 
particle  measurements.  For  example,  Selesnick  et  al.  [1997] 
have  found  that  within  the  interval  of  L  =  3  to  6,  the 
observed  dynamics  of  3  to  8  MeV  electrons  over  a  3-month 
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period  was  consistent  with  a  DLL  =  1.9  •  10-I°  Lll'7±1'3. 
This  compares  to  our  result  for  p  =  5000  MeV/G,  Kp  =  3, 
which  may  be  fit  by  the  form  Dll  =  9.7  •  10-10  L10,5. 
At  L  =  4.5  (corresponding  to  an  energy  of  3.7  MeV)  the 
DLl  from  Selesnick  et  al.  evaluates  to  8.3  •  10-3  day-1, 
which  is  comparable  to  our  value  of  7.0  •  10-3  day-1. 

6.  Summary 

[45]  Using  the  CRRES  electric  field  data  we  have  modeled 
the  power  spectral  density  P(f)  for  frequencies  between  0.2 
and  16  mHz,  as  a  function  of  L  and  Kp;  the  dependence  on 
local  time  was  not  examined.  A  significant  proportion  of  the 
spectra  exhibit  strong  peaks  above  2  mHz  (most  predomi¬ 
nantly  at  ~4  mHz),  thus  precluding  a  meaningful  power  law 
fit  at  higher  Kp  as  was  performed  by  Mozer  [1971]  and 
Holzworth  and  Mozer  [1979].  Instead,  the  spectral  power 
was  fit  to  the  function  P(L,  Kp)  =  a  Lb  exp(cKp)  at  each 
frequency  independently. 

[46]  The  model  spectra  (in  log  power)  are  characterized 
by  a  standard  deviation  of  ~1  (Figures  3  and  4)  meaning 
that  more  than  16%  of  the  individual  power  spectra  differ 
from  the  model  power  spectra  by  a  factor  of  ~  10  as  seen  in 
Figure  8.  This  also  means  that  the  true  characterization  of 
the  environment  may  require  Dll  to  be  at  least  a  factor  of 
10  higher  or  lower  than  the  modeled  Dll  despite  the  fact 
that  they  have  been  parameterized  by  Kp  to  incorporate  the 
short-term  variability. 

[47]  There  are  several  significant  uncertainties  involved 
with  the  assumptions  necessitated  by  single  point  measure¬ 
ments.  Lanzerotti  and  Wolfe  [1980]  concluded  their  study 
comparing  empirical  estimates  of  Dll  and  DlL  by  stating 
that,  “Substantial  more  work  in  comparing,  for  the  same 
time  intervals,  electric  and  magnetic  field  fluctuations  in 
space,  on  the  ground,  and  at  different  latitudes  is  required  in 
order  to  ultimately  completely  understand  the  temporal  and 
spatial  dependencies  of  the  third-invariant-violating  pro¬ 
cesses  in  the  earth’s  magnetosphere.”  Nearly  a  quarter  of 
a  century  later,  we  find  ourselves  in  much  the  same 
situation.  This  study  has  taken  a  step  forward  by  estimating 
the  Dll  based  on  single  point  measurements.  However, 
given  that  it  is  the  global  structure  and  variability  of  the 
electric  and  magnetic  fields  that  control  the  radiation  belt 
particle  populations,  it  will  remain  near  impossible  to 
disentangle  the  various  competing  acceleration  and  loss 
mechanisms  until  we  can  better  define  what  those  fields 
are  doing  through  the  deployment  of  multisatellite  constel¬ 
lations  in  coordination  with  ground  stations. 
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